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Abstract. 

We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in 
comparison to previous linear perturbation analysis of Kelvin-Helmholtz insta- 
bility developing in the jet. Our aim is to assess advantages and limitations of 
both analytical and numerical approaches and to identify spatial and temporal 
scales on which the linear regime of Kelvin-Helmholtz instability can be applied 
in studies of morphology and kinematics of parsec-scale jets. 
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1. Introduction 

3C 273 is the brightest quasar known. Due 
to its relative closeness (z — 0.158), it 
has become a prime target for studies in 
all spectral bands aimed at understand- 
ing t he AGN phenomenon ijCourvoisier I 
1998). Radio emission from the jet in 
3C 273 has been observed on parsec 
scales using VLB ! 1 llPearson et al. M198H 
Krichbaum et al ~1 1200(1 : lAbraham et al. I 
Tflflfit iLohanov fc Zen^T 2001). 
Lobanov fc Zensus I l)2001l LZOl here- 
after) interpreted structures observed in 
the parsec-scale jet of 3C 273 as a double 
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helix and applied linear ana lysis of Kelvin- 
Helmholtz instability (e.g., lHardeel l2000|) 
to infer the basic p arameters of the flow. 
lAsada et al. I l|2002t suggest that a helical 
magnetic field could generate similar 
patterns in the jet. LZOl obtain a bulk 
Lorentz factor W = 2.1, well below those 
measured for superluminal compo nents 
(W = 5 - 10, lAbraham et~aT~l ll996Jb but 
suggest that instabilities develop in the 
slower underlying flow, interpreting those 
fast components as shock waves inside the 
jet. This result, if confirmed, would make 
the combination of observations and linear 
stability theory a powerful tool for probing 
the physical conditions in relativistic jets. 
We aim here to compare the results of 
the analytical modelling to the structures 
generated in a numerical simulation of a 
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steady jet, with the initial conditions of 
the underlying flow given in LZ01. The 
jet is perturbed with several helical and 
elliptical modes, which correspond to the 
modes identified by LZ01. Previous works 
by Perucho et al. l)2004at l2004b|) have 
shown that numerical simulations can be 
used to study the transition from the linear 
to the non-linear regime, and we will apply 
this approach to investigate whether the 
results obtained from the numerical and 
linear methods can be reconciled. 

2. Numerical simulation 

2.1. Initial setup 

We start with a steady jet, which has 
a Lorentz factor W — 2.1, density con- 
trast with the external medium r\ = 0.023, 
sound speed c s _j — 0.53 c in the jet and 
c s .a — 0.08 c in the external medium and 
perfect gas equation of state (with adia- 
batic exponent 7 = 4/3). Assuming an 
angle to the line of sight 8 = 15° and 
redshift z = 0.158 (lmas = 2.43 pc), 
the observed jet is 169/i _1 pc long (with 
the modified Hubble constant Ho — 70 x 
ft,kms _1 Mpc _1 . Considering the jet radius 
given in LZ01 (0.8 pc), the numerical grid 
is 211 Rj (axial) times 8 Rj (transversal), 
i.e., 169 pc x 6.4 pc. 

Resolution is 16 cells/_Rj in the 
transversal direction and 4cells/i?j in the 
direction of the flow. A shear layer of 
2 Rj width is included in the initial rest 
mass density and axial velocity profiles to 
keep numerical stability of the initial jet. 
Elliptical and helical modes are induced at 
the inlet. 

Frequencies of the excited modes are 
computed from the observed wavelengths, 
\ obs , corrected for projection effects, rel- 
ativistic motion and propagation speed 
of the instability, v w . This yields oj — 
2irv w /\ theor , with 

theor _ X obs (l-v w /c cosg) 
A - skrtf ' (1) 

We use v w = 0.23 c for all body modes 
identified in LZ01 (see Table 1). The sur- 



face modes are assumed to mov e with a 
speed close to that of the jet (e.g.. lHardee I 
2000). The 18 mas mode in LZ01, which is 
not identified with any Kelvin-Hclmholtz 
mode, was not included in this simula- 
tion. Table 1 lists the modes which have 
been excited in the simulation. The sim- 
ulation lasted for a time 1097 Rj/c (i.e., 
~ 2852 yrs), and it used ~ 11 Gb RAM 
Memory and 8 processors during around 30 
days in a SGI Altix 3000 computer. 

2.2. Discussion 

The simulation run can be divided in two 
different parts. In the first part, the insta- 
bility grows linearly. In the second part, 
disruption occurs and it dominates further 
evolution of the flow. In Figs. and we 
display several transversal cuts and axial 
cuts at two different times of the simula- 
tion, respectively. The transverse cuts illus- 
trate mode competition showing that dif- 
ferent modes dominate at different posi- 
tions and times in the jet. The longitudinal 
structures presented in Figs. El and indi- 
cate that the linear phase of the instabil- 
ity growth is dominated by the short heli- 
cal mode (~ 4 Rj) modulated by the longer 
helical (antisymmetric, ~ 20 Rj) and ellip- 
tical (symmetric, ~ 50 Rj) modes. 

The linear growth of instability ends up 
with the disruption of the jet caused by one 
of the longest helical modes (~ 20 Rj), at 
time t = 350 Rj/c. After disruption, evolu- 
tion of the jet is totally influenced by this 
mode, and it induces perturbations that 
propagate slowly backwards as a backflow 
in the ambient medium. The disruption 
point moves outwards due to the constant 
injection of momentum at the inlet and the 
change of conditions around the jet, which 
seem to make it more stable. This point 
moves from 160 Rj to 180 Rj by the end of 
the simulation (see Fig. EI> • 

To check the consistency of the results 
obtained for the linear regime of the simu- 
lation with observed structures, we need to 
measure the propagation speeds of the per- 
turbations. The time step between subse- 
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A obs 
(mas) 


Mode 


m 




(mas) 


(mas) 


(mas) 
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Hb2 


18.7 


4 


0.44 


0.54 


2.27 


4 


Ebi , Hbi 


37.4 


25 


2.7 


3.37 


14.3 


12 


E a 


21. T 


50 


5.5 


6.7 


28.5 



Table 1. First two columns give identified wavelengths and modes in LZ01 (H stands for 
helical, E for elliptical, s for surface mode and 61 and 62 for first and second body modes, 
respectively), third column gives the intrinsic wavelengths (see text), in the fourth column 
we have written observed wavelengths in the simulation, and the last three columns give 
the fourth column wavelengths as observed depending on the wave speed. a Computed 
assuming it propagates with the flow speed. 
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Fig. 1. Map of pressure perturbation transversal cuts in arbitrary units, values increasing 
from dark to lighter colors. Solid line indicates v z = 0.8 c contour. Three left panels: Cuts 
at 35i?j, t — 70, 140, 200 Rj/c where elliptical mode rotation is apparent. Three right 
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Fig. 2. Map of Lorentz factor distribution of a portion of the jet at a time before disrup- 
tion, where a large amplitude wave is apparent (top panel, t = 320 Rj/c) and at the last 
frame (bottom panel, t — 1097 Rj/c). Coordinates are in jet radii. Note that the vertical 
scale size is increased by a factor of 4 for the sake of clarity. 



quent frames is too large to measure these 
speeds directly from the simulation data, 
and we use pressure perturbations profiles 
for this purpose. From pressure perturba- 
tion plots (see Fig.|3J), we estimate the ve- 
locity of propagation of the fastest modes 
by measuring the position of the perturba- 
tion front in each frame. In this way, we 
find perturbations which travel with a ve- 



locity close to that of the flow (v w ~ 0.88 c 
as an upper limit). We can also derive 
the wave speed of the disruptive mode fol- 
lowing the motion of the large amplitude 
wave (see Fig. [5| from frame to frame, 
and we find v w ~ 0.38 c. We can asso- 
ciate the former, faster perturbation, with 
the longer wavelength and longer exponen- 
tial growth length elliptical surface mode 
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Fig. 3. Longitudinal cut of pressure per- 
turbation at Rj/2 in symmetric positions 
with respect to the jet axis at t = 250 Rj/c. 
Helical (antisymmetric) structures of 4 and 
25 Rj and an elliptic (symmetric) one of 
50 Rj are apparent. 



(see Fig. |3J), and the latter with a shorter 
wavelength and shorte r exponential g rowth 
length body mode fsee lHardee 1200(1 ). Both 
measured speeds are different to that given 
by LZ01 (0.23 c). This difference may be 
caused by an accumulation of errors in dif- 
ferent assumptions, or even by the excita- 
tion of a slightly different (implying dif- 
ferent velocity). The observer frame wave- 
lengths calculated for different values of the 
propagation speed are given in Table 1. 

3. Conclusions 

We have shown that the structures found in 
a jet with the physical properties of the un- 
derlying flow given in LZ01, and perturbed 
with elliptical and helical modes, are of the 
same order in size as those observed, if rela- 
tivistic propagation effects of the waves are 
taken into account. Although there is a dif- 
ference between wave speeds found in this 
work and that derived from linear analy- 
sis and we do not find a unique combina- 
tion of parameters which explain observed 
structures, we show in Table 1 that simi- 
lar wavelengths to those observed are found 
for given combinations of measured wave- 
lengths in the simulation and wave speeds. 



Our simulations do not address the 
question of long term stability of the jet. 
The simulated jet is disrupted at distances 
of ~ 170 hT 1 pc, which is much shorter than 
the observed ~ 60/i -1 kpc extent of the 
jet in 3C273. A combination of five differ- 
ent factors may be responsible for this dis- 
crepancy. 1) Magnetic fields have not been 
taken into account neither in the linear 
analysis, nor in the numerical simulation, 
but they may be dynamically important at 
parsec scales. 2) We only simulate the un- 
derlying flow, without taking into account 
the superluminal components. 3) Initial 
amplitudes of perturbations cannot be con- 
strained accurately. 4) Inaccuracies in the 
approximations used in the linear analysis 
can lead to differences in physical parame- 
ters derived. 5) Differential rotation of the 
jet, shear layer thickness, and a decreasing 
density of the external medium could affect 
the wavelengths and speeds of the insta- 
bility modes. The combined effect of these 
factors could well change the stability prop- 
erties of the jet. 
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